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ABSTRACT 



The electrical resistivity of the accreted mountain in a millisecond pulsar is lim- 
ited by the observed spin-down rate of binary radio millisecond pulsars (BRM- 
SPs) and the spins and X-ray fluxes of accreting millisecond pulsars (AMSPs). 
We find 7] > 10~ 28 s (tsd/1 Gyr) -0 ' 8 (where tsd is the spin-down age) for BRMSPs 
and rj > 10~ 25 s (M^/Me) ' 6 (where M a and Me are the actual and Eddington ac- 
cretion rates) for AMSPs. These limits are inferred assuming that the mountain 
attains a steady state, where matter diffuses resistively across magnetic flux sur- 
faces but is replenished at an equal rate by infalling material. The mountain then 
relaxes further resistively after accretion ceases. The BRMSP spin-down limit 
approaches the theoretical electron-impurity resistivity at temperatures > 10 5 K 
for an impurity concentration of ~ 0.1, while the AMSP stalling limit falls two 
orders of magnitude below the theoretical electron-phonon resistivity for temper- 
atures above 10 8 K. Hence BRMSP observations are already challenging theoret- 
ical resistivity calculations in a useful way. Next-generation gravitational-wave 
interferometers will constrain t] at a level that will be competitive with electro- 
magnetic observations. 

Subject headings: accretion, accretion disks - stars: magnetic fields - stars: neutron - 
pulsars: general. 
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Introduction 



R ecent searches for periodic gravitation al wave signals f rom iso 



20071); 



Wette et al 



2008 



Abbott et al. 



20081 ) and accreting ((Abbott 



ated (Abbott et al 



2007b 



Watts et al 



20081 ) neutron stars with the Laser Interferometer Gravitational Wave Observatory (LIGO) 
have placed new upper limits on the mass quadrupole moments of these objects. In rare 
instances, the LIGO upper limit beats the indirect limit inferred from electromagnetic 



spin-down observations, e.g. the Crab's mass ellipticity is measured 
e < 1.8 x 10~ 4 , approximately 25% of the indirect limit (jAbbott et al 



)v L 



GO to be 



20081 ). More often, 



the indirect limit bea ts the LIGO limit, e.g. for slowly decelerating radio millisecond pulsars 



flAbbott et al. 



2007bJ). 



As theoretical models for the quadrupole become more accurate and sophisticated, 
spin-down measurements and LIGO nondetections translate into stricter limits on the 
constitutive properties of a neutron star. For example, the quadrupole produced by 
electron-capture gradients requires a persistent lateral temperature difference of ~ 5% at 
the base of the outer crust, and a normal rather than superfluid core, to be consistent 
with the gravitational-w ave stalling interpretation of low-mass X-ray binary (LMXB) spins 



rtUshomirsky et al. 



20001 ) . The compressibility of nuclear matter must be lower than ~ 200 



MeV to explain the nondetec tion of a gravitational wave signal from large pulsar glitches 



( Ivan Eysden fc Melatos 



20081 ). The breaking strain of the crust is bounde d above by th e 



nondetection of a gravi t ational- wave signal from accreting neutron stars (IHaskell et al. 



2007 



Watts et al 



2008 



Horowitz fc Kadau 



2009|) 



In this paper, we show how the observed spin-down rates of binary radio millisecond 
pulsars (BRMSPs), and the spins and X-ray fluxes of accreting millisecond pulsars (AMSPs) 
translate into lower bounds on the electrical resistivity of an accreted neutron star crust. 
The bounds arise because the magnetically confined polar mountain formed during 
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accretion (jPayne fc Melatosl 12004 ; iMelatos fc Paynd 120051 ) must relax resistively, until the 
quadrupole is low enough to respect the observational limits. Such an argument can be 
mounted now with rising confidence, because there have been substantial improvements 
recently in our understanding of the structure of magnetic mo untains, including 



studies of ideal-ma 



gneto 



Vigelius fc Melatos 



equati ons of state 



crust (jWette et ai- 



r y dro dy namic (ideal-MHD ) instabilities ([Pavne fc Melatos 



2008 



2009af). Ohmic diff usion ( jVigelius fc Melatos 



2006 



Privmak fc Melatos 



2009^), realistic 



20091 ). and hydrodynamic sinking into the non-rigid 



20091 ). The theory is now robust enough that its prediction of large 
ideal-MHD quadrupoles (e ~ 10~ 5 ) deserves to be taken seriously. One must therefore 
appeal to nonideal processes, like resistive relaxation, to bring e below ~ 10~ 8 , as the 
observational limits require. 

The paper is set out as follows. In §2J magnetic mountain models are briefly reviewed, 
and an approximate formula is given for e as a function of accreted mass, electrical 
resistivity, and age. The formula is applied to convert observations of BRMSP spin-down 
rates (§3J) and AMSP spin frequencies (§4|) into lower limits on the electrical resistivity, 
assuming energy conservation and gravitational- wave stalling respectively. In §5j the 
limits are mapped onto a resistivity-versus-temperature plot and compared with the 



latest theoretical calcu 



( lltoh fc Kohyama 



1993 ; 



ations of the electron-phonon and electron-impurity resistivities 



Baiko fc Yakovlev 



19961 



Cumming et al. 



2001 



20041 ). The prospect 



of better limits from future LIGO experiments is canvassed briefly. 



2. Magnetic mountain quadrupole 

The quadrupole generated by hydromagnetic stresses from the natal magnetic field of 
a neutron star (before accretion begins) is too small to emit gravitational waves at a level 
detectable in the spin down of MSPs or by LIGO. For dipolar fields in the range 10 12 G 
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< B < 10 13 G (jHartman et al.lll997t lArzoumanian et al.ll2002l ; iFaucher-Giguere fc Kaspi 
20061 ). one obtains ho < 10 -30 for the gravitational wave strain and |/*| <4x 10 -32 Hz s" 1 



for the frequency spin-down rate (IBonazzola Sz Gourgoulhonlll996l ). A stochastic magnetic 
field in a type I supercon ducting core , or a strong toroidal field (> 10 14 G), give h ~ 10~ 27 



and |/„| ~ 10~ 25 Hz s^ 1 ([Cutler 



2002) 



On the other hand, magnetic burial by accretio n leads to a much 



l arger 



Payne fc Melatosl (12004 ) computed 



quadrupole, 



as accreted material is funneled onto the polar cap. 
the unique, quasistatic sequence of ideal-MHD equilibria that describes how burial 
proceeds as a function of the total accreted mass M a while self-consistently respecting the 
flux-freezing constraint! 1 ]. They found that the accreted matter is confined at the poles by 



an equatorial belt of compressed magnetic field a nchored in the deep crust (jHameury et al. 



1983 



Brown fc Bildsten 



1998 



Litwin et al 



20011 ). The resulting 'magnetic mountain' 



has a substantial quadrupole moment, even after it relaxes hydro magnetically through 
the undular Parker inst ability. Expressed as a mass ellipticity eMHn, the final ideal-MHD 



quadrupole is given by ( 



2006 



Vigelius fc Melatos 



avne fc Melatos 



2009a; 



2004; 



Wette et al. 



Melatos fc Payne 



20091 ) 



2005 



Payne fc Melatos 



£mhd 



Ma ( 1 Ma 

M r \ M r 



2 x 10- 7 

with M c » 2 x 1O~ 5 M . Numerical simulations, where a mountain is grown ab initio by 
mass injection, confirm the overall scaling for eMHD [with some extra oscillatory dynamics 



f lVigelius fc Melatosl l2009a{ )]. 



Magnetic mountains relax resistively on the diffusive time scale, Td = 8.6 x 10 7 yr, 



1 Flux- freezing is valid in the crust, down to a depth where carbon ignition occurs in 
one-dimensional models, but not necessarily in the magnetic belt region, where the field is 
highly distorted and diffusion proceeds faster. 
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after accretion switches off (IVigelius Sz Melatos! l2009bf ). In a BRMSP, resistive relaxation 
proceeds roughly exponentially, with 



e = eMHDexp(-t?7/r d ) 



(2) 



where r\ measures the resistivity r\ in uni ts of 1.3 x 10 7 s, and t is the time elapsed 



since accretion stops [compare Fig. 1 in 



Vigelius & Melatosl (J2009bJ)]. In an AMSP 



where accretion continues today, a steady state is established, where the influx of accreted 
material into a magnetic flux tube exactly replenishes the efflux due to Ohmic diffusion. 
The steady state is attained when the diffusion time-scale equals the accretion time-scale. 
As 7~ri is directly proportional to 77, the correspond ing saturation ellipticity scales as 



f lBrown Sz Bildsten 



1998 ; 



Vigelius Sz Melatos 



2009af ) 



e = min e MH D, 5.1 x lO'V M a /M E 



where M a and Me are the actual and Eddington accretion rates, respectively. 



(3) 



Mountain relaxation is required by existing data. LMXBs s hould have been detected b y 



LIGO already, if the quadrupole is as large as Eq. 



Yet directed searches have found nothing to date (jAbbott 



([I]) predicts fjyigelius Sz Me 



2007b 



Abbott et al. 



atos 



2009ah . 



200731 ). 



In general, the magnetic mountain axis is tilted with respect to the rotation axis. 
Hence the star precessesj, generat ing gravitational waves at f* and 2/*. The polarized way e 



strain is written down elsewhere ( 



Van Den Broeck 



2005 



Zimmermann Sz Szedenits 



Vigelius Sz Melatos 



1979 



Jaranowski et al. 



1998 



2009al ). Its characteristic amplitude is 



c 4 V 



(4) 



2 It is therefore possible in principle to measure independently the quadrup ole moment 



from modulations in X-ray light curves, e.g. in RCW 103 or XTE J1814— 338 fjHeyl 



2002 



Chung et al. 



20081). 
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where I zz is the moment of inertia, /* is the spin frequency, D is the distance to the source, 
9 is the wobble angle, and % is the inclination angle of the line of sight. For small 6, we have 
F(6,i)<x6. 

We combine €mhd from Eq. ([1]) with the resistive decay prescriptions in this section 
and Eq. (J3J) to derive limits on 77 in BRMSPs and AMSPs in §3] and §H] respectively. Ideally, 
e for both BRMSPs and AMS Ps should be estimated by the same method , since every 



AMSP evolves into a BMRSP f lWiinands fc van der Klis 



1998 ; 



Alpar 



2008 



But we are 



prevented from making this connection because two crucial pieces of important information 
are missing: M a for BRMSPs, and the age of AMSPs. 

3. Binary radio millisecond pulsars 

3.1. Spin down 

If we attribute the loss of angular momentum observed in pulsar timing experiments 
to gravitational wave emissi on alone, we can wr i te down a corresp onding limit on the 



gravitational wave strength ( lAbbott et al. 



2007b 



Wette et al. 



2008|) 



hsD < 



5GI zz \f\ 
2c 3 D 2 f 



1/2 



(5) 



which translates to an ellipticity 

e SD < 5 x 10~ 8 



l/l 



1/2 



lO-^Hzs- 1 



/ 



100 Hz 



-5/2 



(6) 



where / is the gravitational wave frequency (assumed here to be 2/*), and / is its time 
derivative. 

We compare cmhd and esD for those BRMSPs in the literature, whose M a has been 
estimated with the help of binary evolutionary simulations assuming that mass transfer is 



( 



Object /* |/,| D M a 
[Hz] 10~ 16 Hz s- 1 [kpc] [M Q ] 

B0021-72E 283 

B0655+64 5.11 

(B0820+02) 1.16 

(B1800-27) 2.99 

(B1831-00) 1.92 

B1913+16 16.0 

B1953+29 163 

J0034-0534 533 

J0218+4232 430 

J0437-4715 174 

J1022+1001 60.0 

J1045-4509 134 

J1713+0747 219 

J2019+2425 254 

J2145-0750 62.0 

J2229+2643 336 

J2317+1439 290 

able 1: Binary radio millise cond pulsars for which M a has been estimated 
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Eddington-limited and conservative ( Ivan den Heuvel fc Bitzarakilll995l ). The objects in this 
class are listed in Table [TJ together with their spin-down parameters, distances, a nd M a 



values. Typical fractional uncertainties in the quoted M a are < 50% flWiiers 



becau se it is unclear how conservative mass transfer really is ( 



1997 ), mainly 



van den Heuvel fc Bitzaraki 



19951 ). The spin parameters and distance, inferred from the dispe rsion measure, are drawn 
from the Australia Telescope National Facility pulsar database! (IHobbs et al.ll2004l ). We 
find that the ideal-MHD ellipticity is higher than the spin-down limit for all but three 
objects (PSR B0820+02, PSR B1800-27, and PSR B1831-00, in parentheses in Tabled]). 
The latter objects are excluded from the subsequent analysis because they do not translate 
into limits on 77. 

Assuming that resistive relaxation reduces eMHD below eso, and that the time available 
to do so is the spin-down age at the present epoch, tsd = — /*/(2/*), we can combine ([I]), 
([2]), and ([6]) to get a minimum resistivity r/ min . The results are presented in Fig. [TJ which 
displays ?7 m i n as a function of tsd on a log-log plot for the objects discussed in the previous 



paragraph (filled squares). A least-squares fit yields r^ 



10 



-27.3±0.8 



tsd/1 Gyr) 



-0.7±0.1 



where we omit the single obvious outlier at tsd = 4.5 x 10 9 yr from the fit. 



20091 ) . If we assume that the 



Equation ([6]) is almost certainly too conservative (lOwenl 
electromagnetic torque exceeds the gravitational-wave torque in BRMSPs by at least the 
same factor ( ~ 4.1 times) as in the Crab pulsar, as indicated by the results of the LIGO 



Crab search ( Abbott et al. 



20081 ) (an arbitrary assumption, to be sure, but plausible), then 



we arrive at a more stringent limit on rj. Replacing \f*\ by 4|/„,| in Eq. (J5J), we find that 
esD doubles for all the objects plotted in Fig. [1] (open circles). Fitting the open circles, we 
obtain Vmm = 10- 27 - 2±a6 s(r sd /l Gyr)-°- 8±0 - 2 



5 www. atnf.csiro.au/research/pulsar/psrcat/ 
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Fig. 1. — Electrical resistivity rj (in s), inferred from the spin-down limit for BRMSPs, 
plotted as a function of characteristic age tr d (in yr). Lower limits on t] from the 14 objects 



listed in Ivan den Heuvel fc Bitzarakil (119951 ) are shown as filled squares; the three objects 
with 6mhd < £sd do not yield useful limits. A more stringent limit on rj (open circles) can be 
obtained by assuming that the electromagnetic torque exceeds the gravitational-wave torque 
by at least a factor of four, like in the Crab (see text). The solid and dashed lines are 
power-law fits to the filled squares and the open circles, respectively. 
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3.2. Nondetection of gravitational waves 

In a similar fashion, we can obtain a less stringent lower limit on 77 from gravitational- 
wave observations. Directed searches for gravitational waves from kno wn radio pulsar s 



set an upper 



Abbott et al. 



limit o n the wave strain (and hence e) for these objects ((Abbott 



2007b ; 



2007al ). We then apply the same argument as in §3.1 



Fig. [2] displays the expected gravitational wave strain h as a function of gravitational- 
wave frequency in the ideal-MHD limit (open diamonds), calculated from eMHD, along 
with the detection thresholds for Initial and Advanced LIGO (dashed and dotted curves 
respectively), where we assume 14 days of coherent integration, a false alarm rate of 1 per 
cent, and a false dismissal rate of 10 per cent. For comparison, we plot /isd (filled squares, 
where the arrows indicate that it is an upper limit). In the ideal- MHD limit, nine obje cts 



2007tJ) 



lie above the Initial LIGO threshold yet have not been detected (lAbbott et al. 
confirming that the ideal-MHD mountain relaxes substantially after accretion shuts off. 
More significantly, for all nine objects, hsn lies below the Advanced LIGO noise curve, 
although PSR J0437— 4715 is close. Next-generation (e.g. subterranean) interferometers 
will significantly lower the detection threshold, to the point where future gravitational-wave 
observations will constrain ?7 min more tightly than spin-down data (which are likely to be 
dominated by the electromagnetic torque). 



4. Accreting millisecond pulsars 

For AMSPs, we can derive a limit on the electrical resistivity by assuming that all 
objects are in spin equilibrium, such that the gr avitational radiatio n react ion torque exactly 



balances the accretion torque (IBildsten 



19981 ) . 



Watts fc Krishnanl (120091 ) cautioned that 



spin balance may not hold on short time-scales. However, even if the torques do not balance 
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Fig. 2. — Expected gravitational w ave strain ho in the ideal-MHD limit for the objects from 



van den Heuvel fc Bitzarakil (119951 ) (open diamonds) and the spin-down lower limit hsn from 
Eq. (jHJ) (filled squares). We also plot the detection threshold, assuming 14 days of coherent 
integration, a false alarm rate of 1 per cent, and a false dismissal rate of 10 per cent, for 
Initial LIGO (dashed) and Advanced LIGO (dotted). Note that, although e saturates for 
all BRMSPs (M a 3> M c ), the expected gravitational wave strain depends on the distance 
(ho oc D~ l ) and is therefore different for individual objects. 
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Fig. 3. — Electrical resistivity r\ (in s), inferred from the stalling frequency of AMSPs, plotted 
as a function of accretion rate M a (in units of the Ed dington rate Me ), inferred from the 



Galloway! (12008 



are shown as filled 



X-ray flux Fx- Lower limits on rj for 17 objects from 
squares (persistent sources) and filled diamonds (transient sources). The open diamonds are 
lower limits inferred by the BRMSP method in §3.11 from the observed spin down between 
outbursts of XTE J0292-314 and SAX J1808.4-3658. The solid (dashed) line represents a 
power-law fit to the persistent (transient) objects. 
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instantaneously, it is enough for the present application that they balance on average over 
the long term. 

Torque balance implies an indirect upper limit on the gravitational wave strain as a 
function of X-ray flux Fx and spin frequency, 



< 



\\ 



1/2 



4x 10~ 27 "" \10 _8 ergs~ 1 cnr 2 
or, equivalently, on the gravitational mass ellipticity 

F X \ 1/2 ( 



U 



300 Hz 



-1/2 



(7) 



< 



/. 



-5/2 



2 x io~ 7 ~ v 10 ~ 8er g s " lcm ~ 2 y V lk P c 7 V30oh Z/ 

if the star is located at a distance D. When interpreting Eq. (jHJ) as an upper limit, it is 
important to keep two things in mind, (i) The X-ray emission is generally anisotropic, (ii) 
The accretion torque includes a magnetic component, whos e lever arm, th e Alfven radius, 



can extend several stellar radii fjMelatos fc Owen 



2009|) [cf. iBildstenl f|l998h ]. 



We are now in a position to compare the resistively relaxed quadrupole given by Eqs. 
(fT)) and ( 3 ) to t he observational upper limit (E]) and thus obtain a lower bound on 77. 
Gallowayl ( 120081 ) tabulated D and Fx for 17 AMSPs. The key properties of these objects 
are listed in Table El Four of the objects are persistent X-ray sources; the remainder are 
transient, with a range of outburst and recurrence times. Burst sources are marked by 
an asterisk. We compute 77 min for these objects from Eqs. ([1]), fl3]), and ([8]), and plot it 
against M a in Fig. [3j where we assume that all mechanical energy of the infalling matter 
is converted into X-ray luminosity, viz. Fx ~ (GMM a /R*)/4:7rD 2 . A power law is a 
reasonable fit to the data. We find r/ min = 10~ 26 - 9±a5 (M a /M E ) 0,9±a2 for the persistent 
objects (solid curve, fitted to the filled squares) and 77 min = 10~ 28±4 (M a /M E )°' 5±0 ' 2 for the 
transient objects (dashed curve, fitted to the filled diamonds). 



wo AMSPs are observed to spin down secularly between outbursts. iGallowav et al. 



( 120021 ) reported /* = (-9.2 ± 0.4) x 10" 14 Hz s" 1 in XTE J0929-314, while 



Hartman et al 
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( 120081 ). who followed the spin evolution of SAX J1808.4— 3658 over seven years and four 
transient outbursts, measured /* = —(5.6 ± 2.0) x 10~ 16 Hz s" 1 . If we conservatively 
attribute the spin down to gravitational wave emission alone, we can obtain an independent 
upper limit on e (and hence ?7 m i n ) by the BRMSP method in §3.11 For comparison, we 
include these lower bounds as open diamonds in Fig. |3j For XTE J0292— 314, the spin-down 
limit falls two orders of magnitude short of the limit provided by torque balance. However, 
in the case of SAX J1808.4— 3658, spin down during quiescence provides a slightly more 
stringent limit than torque balance. 



5. Resistivity versus temperature 



In this section, we compare the observational limits derived above with theoretical 
estimates for the crustal resistivity. As t] is sensitive to the temperature T, which is different 
in BRMSPs and AMSPs, the comparison covers two decades of temperature and hence a 
range of scattering mechanisms. 

The thermal relaxation time of a neutron star (~ 10 4 yr) is typically shorter than 
the accretion ti me-scale (~ 10 6 yr) . Hen ce, the star is in thermal equilibrium during its 



X-ray lifetime f lYakovlev &: Pethick 



20041 ). with its internal temperature set by M a . Once 



accretion stops, and a BRMSP forms, the crust cools through a combination of neutrino 
(e.g. via direct Urea processes) and ph oton cooling. Competing mechan isms, such as 



viscous diss ipation of rotationa l 



dissipation ( 



Haensel et al. 



heating ( IReisenegger 



1995 



1990 



energy (IPage 



1998 



Yakovlev et al. 



Cheng fc Dai 



1996 



1999 



Yakovlev et al. 



Miralles et al. 



Reisenegger 



1997 : 



1998 



1999), magnetic field 



) and rotochemical 



Iida fc Sato 



19971 ) act to 



reheat the crust. In this paper, for definiteness, we assume that reheating occurs mainly 
rotochemically: the rate of reactions restoring chemical equilibrium is slower than the rate 
of change of particle concentrations during spin down (as the centrifugal force diminishes), 
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Object 



D 

[kpc] 



/. 



M a . 



nu ,™-n 



[Hz] [10- n M o yr 



Persistent sources: 








4U 1916-05* 


8.9 


270 


22 


4U 1702-429* 


5.5 


329 


35 


4U 1728-34* 


5.2 


363 


45 


Transient sources: 








MXB 1658-298* 


12 


567 


60 


Aql X-l 


5.0 


550 


110 


XTE J1814-338 


8.0 


314 


1.4 


4U 1608-52* 


4.1 


620 


100 


XTE J0929-314 


4.0 


185 


0.40 


Swift J1756.9-2508 


8.0 


182 


0.14 


XTE J1807-294 


5.0 


191 


0.30 


GRS 1741.9-2853* 


8.0 


589 


0.010 


XTE J 175 1-305 


8.0 


435 


1.2 


IGR 00291+5934 


5.0 


599 


1.3 


4U 1636-536 


6.0 


581 


55 


HETE J1900. 1-2455 


5.0 


377 


7.0 


SAX J1808.4-3658 


3.4 


401 


0.60 


SAX J1748.9-2021 


8.1 


442 


8.0 


EXO 0748-676* 


7.5 


45.0 


20 


SAX J1750.8-2900* 


6.8 


601 


20 


A 1744-361* 


9.0 


530 


230 


KS 1731-26* 


7.2 


524 


85 



Table 2: Known accreting millisecond rotators [adapted from iGallowayi (120081 ) and 



Watts et al. 



(120081 )] . Type I burst oscillation sources are marked with an asterisk. 
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so that the star is permanently out of chemical equilibrium and generates heat. If the 
spin-down timescale exceeds other time-scales, the temperatu re depends only on the current 



spin-down parameters, with ( Fernandez Sz Reiseneggei 

2/7 



2005|) 



r 

T 



I/. 



10- 16 Hzs- 



/* 



100 Hz 



2/7 



(9) 



and 2.8 < T n /10 4 K < 4.2. The range in T n encompass e s two realistic equa tions of state for 



the crust (IHaensel fc Pichon 



of state for the core (jPrakash et al. 



1994 



Pethick et al 



whole star (jShapiro fc Teukolsky 



1988 



1995 ; 



Akmal et al. 



19981 ) . five equations 



a nd a non-interacting Ferm i -gas m odel for the 



19831 ); see Fernandez &: Reiseneggerl (120051 ) for details. 



Below, we use Eq. fl9]) to compute T for BRMSPs. 

No simple T-M a relation is available for accreting neutron stars. For the purpose 
of this paper, we are sensitive to t] at the base of the mountain, where the bulk of the 
accreted matter is confined. The base is located ~ 0.1 — 1 km below the surface for 
M a f« 0.1M Q , where we assume that the equation of state in the crust can be mod elled 



using a Skyrme-type effective nucleon-nucleon interaction (IDouchin &: Haensel 



accurate treatment of the sinking problem is presented elsewhere ( jWette et al. 



2001). An 



20091). The 



temperature profile of the accretion column, and hence T at the mountain base, is set by 
nuclear heat sources inside the star, the composition of the infalling matter, and the local 
M a . It is not directly observable. 

We consider two scenarios, (i) For transient sources, the temperature in the outer crust 
is set by a balance between deep-crustal heating during outbursts and subsequent neutrino 



cooling; ( 



Brown et al. 



( iColpi et al 



19981 ) . The core temperature attains a steady state after ~ 10 4 yr 
20001 ) and the time-averaged X-ray flux (F acc ) « 0.2 (M a \ c 2 /(AttD 2 ) can be 



relate d to the quiescent flu x by F q w (F acc ) /135, where we assume an accretion efficiency of 



20% (IRutledge et al 



20021 ). We estimate the temperature by ap plying Fick's law, assum ing 



an average thermal conductivity of k ~ 10 18 erg s 1 cm 1 K 1 ( )Haensel fc Zduniklll990l ) 
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and a thermal length scale of I ~ 1 km. This procedure is only applicable if the duration 



of the outbursts is shorter than the thermal diffusion time of the crust, which is not a 



true. For example, AMSP KS 1731—260 accretes actively for > 13 yr (IRutledge et al. 



Brown fc Cumming 



ways 



2002 



2009|). 



We note that, according to the Wiedemann- Franz law, the thermal conductivity k and 
the electrical conductivity a are related through njo = 7r 2 /3(kB/e) 2 T at a temperature 
T, where e denotes the electron charge in cgs units. This agrees well with the theoretical 
values for a (dashed lines in Fig. H]) under the assumption that n ~ 10 18 erg s _1 cm" 1 K -1 . 

(ii) For persistent sources, most of the heat released by thermonuclear burning 
is immediately radiated outward and the interior thermal balance is dominated by 
pycnonuclear reactions in the deep crust. Thus, the temperature in the photosphere is set 
by the accretion energy flux (~ GMM/r) and gradually increases with depth by a factor 
of 20 — 40 to reach (8 — 70) x 10 8 K at ~ 1 km (for lo cal accretion rates per unit area in 



the range 5 < m/rriEdd < 40) (IBrown fc Bildsten 



19981 1 . In order to place a lower bound on 



T (and hence a lower limit on 77), we adopt two independent approaches. First, we assume 
that the observed X-ray flux is thermal and emitted by a polar cap with area ~ 1 km 2 . 
We then apply the Stefan-Boltzmann blackbody formula, neglecting free-free corrections 



( iZavlin et al. 



19961 1 . The results are verified b y comparing with X -ray spectral models 



20051 ). Second, we follow a 



in several objects, e.g. for IGR J00291+5934 (IFalanga et al. 
procedure similar to the transient sources by assuming that the temperature gradient is 
balanced by the quiescent flux (which is not observed directly but is less than the observed 
flux) and again apply Fick's law. 

We combine the above T estimates with the observational rj limits in Fig. HI which 
displays ?7 m i n (in s) as a function of T (in K). For BRMSPs, our fit for t] m i n as a function 
of tsd in Fig. [H and the T-tsd relation in Eq. fl9]), convert into a resistivity-temperature 
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Fig. 4. — Resistivity 7] (in s) as a function of temperature T (in K). We show the theoretical 
estimates for combined electron-impurity and electron-phonon scattering as solid curves 
(upper curve: impurity concentration Q = 10, lower curve: Q = 0.1) and electron-phonon 
scattering alone (dotted curve). Lower limits r) m i n (T) for BRMSPs at the lower (higher) end 
of the range are given by the range of rotochemical equations of state [Eq. (j9])]. These limits 
are shown as upper (lower) thick curves at the left side of the plot. The dashed curves assume 
that the electromagnetic torque exceeds the gravitational- wave torque by a factor of four. 
Vmm{T) for AMSPs are displayed as thick curves at the right side of the plot. We distinguish 
between transient (left diagonal line) and persistent (right solid and dashed diagonal lines) 
accretors. For the latter, we compute the crustal temperature according to Fick's law (solid 
line) and, independently, by assuming black-body emission (dashed line). The lines are all 
lower limits. 
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relation in the range 2.3 x 10 4 K < T < 1.6 x 10 5 K. We plot rj m i n (T) as a pair of thick solid 
curves in the left half of Fig. HJ the lower and upper curves correspond to T = 4.2 x 10 4 
K and T = 2.8 x 10 4 K respectively [cf. Eq. ©]. For AMSPs, we translate 77 min (M a ) into 
Vmm(T) as detailed above, where we distinguish between (i) persistent accretors assuming 
black body radiation (right-most, dashed diagonal line), persistent accretors applying 
Fick's law (lower diagonal line), and transient accretors (upper diagonal line), in the range 
10 -13 < M a /M Q yr _1 < 10~ 8 populated by the objects in Fig. [31 For comparison, we 
overplot the theoretical resistivity for combined electron-impurity and electron-phonon 
scattering for impurity concentrations Q = 10 and 0.1 (upper and lower dashed curves, 
respectively). 

The electrical conductiv ity of the neutron sta r crust is caused by electron scattering 



off phonons and impurities (jCumming et al. 



2001 



20041 ). For temperatures below the 



Umklapp temperature (Ty ~ 10 7 K), electron-phonon scattering is suppressed. In rapid 
accretors (M a > 10~ ll M Q yr" 1 ), phonon scattering dominates for Q < 1. The impurity 
fraction depends on the composition of the ashes produced in steady-state nuclear burning 



at low densities. 



Schatz et al. 



(119991 ) found a large variety of nu c 



Q ~ 100, except in rapid accretors (Q ~ 1 for M > 30Me). I Jones! (j2004al ) noted that, if the 



ei in the crust, with 



primordial crust is completely replaced by heterogen eous accreted ma t ter, o ne has Q 3> 1 



Other authors adopt lower fractions. For example, 



Pons &: Geppert 



(120071 ) preferred 



10 4 ~ Q ~ 10 2 i R their model for magnetic field dissipation. 



6. Discussion 



Fig. H] indicates that spin-down and indirect gravitational- wave stalling limits are on 
the verge of challenging theoretical models of the electrical resistivity of the accreted crust 
in a neutron star. The limits from BRMSPs assuming rotochemical heating intersect the 
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Q = 0.1 electron-impurity curve above T > 10 5 K. Hence, very young BRMSPs (upper 
end of the curves) tentatively exclude impurity concentrations Q < 0.1 and will place even 
stricter limits on Q as younger objects are found. 



Brown fc Cummingi (120091 ) constrained Q by fitting theoretical cooling models 
to observed post-outburst light curves of KS 1731—260 and MXB 165 9—29. These 



autho rs found Q ~ 1, and certainly Q <C 10, in close agree ment with 



( 2007 ) . Molecular dyna mics simulations confirm this result (IHorowitz et al. 



Shternin et al. 



Horowitz &: Kadau 



2008 : 



20091 ). revealing a phase-separated regular crystal, with low-Q regions 



embedded in a high-Q phase. One interesting possibility is that the solid ph ase of 



the crust 



in isolated neutron stars is heterogenous with respect to the nuclear charge ((Jones 



2004bh . 



lowering the resistivity by four orders of m agnitude to rj ps 10 24 s. Recent ob servations 

Is 



of cooling curves after extended outbursts (jWijnand 



2001; 



Cackett 



et al 



high thermal condu ctivity of the crust, typical of a regular crystal (IRutledge et al. 



Shternin et al 



20070 



20061) suggest a 



2002 



Future gravitational wave experiments will tighten the above limits. Advanced LIGO, 
for example, offers a narrowband co nfiguration, which lowers the instrumental noise by two 



orders of mag nitude (lAbbott 



uncertainties (IWatts et al 



2007aj). The main challenge to searches is posed by parameter 



20081 ). Here, simultaneous electromagnetic observations will 



prove beneficial, reducing the number of search templates and the computational effort. 
In any case, this paper illustrates one of the ways in which even gravitational-wave 
nondetections yield important independent information about the constitutive properties of 
neutron stars. 

Several shortcomings in the theoretical modelling of magnetic mo untains still need 



to be addressed, most notably sinking of the mountain into the crust (jWette et al 



20091 ) 



and a realistic equation of state which includes nuclear reactions. In addition, in the 
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context of this article, it is important to ultimately include the Hall effect. While it is not 



dissipative, the 



iall effect transfers magnetic e nergy from large to small scales through 



a Hall cascade (IGoldreich fc Reisenegger 



subsequently decay on the Ohmic timesca 



for neutron star crusts by 



Cumming et al. 



19921 ) . The small-scale magnetic structures 
e. The Hall timescale, tnaih has been computed 
(bfJOJ), who found 10 3 < t Ha ii/yr < 10 7 , well 



below the typical X-ray lifetime of an LMXB. This result, however, is only valid as long 



as the impurity factor Q is low 



Hall dissipation fjCumming et al. 



or h igher Q > 1, Ohmic dissipation always dominates 



20041 ). Furthermore, a subsurface toroidal magnetic field 
can create small-scale poloidal field structures through a Hall-drift induced instability, 
where a background shear in the electron velocity drives the growth of long-wav elength 



modes (jCumming et al. 



2004 ; 



Geppert et al. 



2003 



Rheinhardt &: Geppert 



20021 ). The 



consequences of this instability for the magnetic field evolution in neutron stars are still 



unclear; perturbations to the crustal field are eff ectively reflectec . 



and might be trapped in the lower crust layers (jCumming et al 



byt 



re solid/fluid interface 



200J). 



In this paper, the theoretical e is intentionally calculated differently for BRMSPs and 
AMSPs. On the face of it, this seems unnecesary, given that BRMSPs evolve from AMSPs. 
But we are forced to take this approach because we do not know the historic M a for any 
specific BRMSP. Hence, we cannot apply Eq. flTJ to get eMHD, while Eq. (J2]) does not apply 
in the context of ongoing accretion. On the plus side, however, the effect of Ohmic diffusion 
during ongoing a ccreti on has not been studied properly in two and three dimensions, as 



Cumming et al. 



(120011 ) pointed out. It is therefore an advantage to have two independent 



calculations for the two classes of objects, pending further studies. 

Finally, if the magnetic mountain picture of field evolution is accurate, then the 
spin-down limits in §3.11 can be combined with magnetic moment data for BRMSPs to test 
the self-consistency of the argument in this paper. For example, if the buried magnetic 
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field is anchored deep in the core, it resurrects as the mountain relaxes resistively, tending 
towards its natal value over many T&. This can lead to an inconsistency if eso/ e MHD is 
too small, for then the resurrected magnetic moment may exceed the value observed. 
Unfortunately, the above test cannot be applied (yet) with confidence, because it rests 
heavily on the assumption that the magnetic field is anchored in the core. If, by contrast, 
the source currents are confined to the crust, minimal field resurrection occurs; both 
e and the magnetic moment decay resistively, and no useful constraint emerges. The 
question of where the field is anchored, crust or core, is hotly debated and remains one 



of the main obstacles in understanding magnetic fie 



( Konar &: Bhattacharya 



1997; 



Ruderman et al. 



1998 



d evolution in acc reting neutron stars 



Pons et al. 



2009J). 
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